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The optimized structure and electronic properties of neutral and singly charged 
magnesium clusters have been investigated using ab initio theoretical methods based 
on density-functional theory and systematic post-Hartree-Fock many-body pertur- 
bation theory accounting for all electrons in the system. We have systematically 
calculated the optimized geometries of neutral and singly charged magnesium clus- 
ters consisting of up to 21 atoms, electronic shell closures, binding energies per atom, 
ionization potentials and the gap between the highest occupied and the lowest unoc- 
cupied molecular orbitals. We have investigated the transition to the hep structure 
and metallic evolution of the magnesium clusters, as well as the stability of linear 
chains and rings of magnesium atoms. The results obtained are compared with the 
available experimental data and the results of other theoretical works. 



I. INTRODUCTION 

Metal clusters have been recognized as new physical objects with their own properties 
almost two decades ago. This became clear after such experimental successes as the discovery 
of electronic shell structure in metal clusters observation of plasmon resonances |2], 
[3], ||, formation of singly and doubly charged negative cluster ions || and many more. 
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Comprehensive survey of the field can be found in review papers and books; see, e.g., 

§ 0, |, |, 0i§§g. 

The electronic shell structure of metal clusters has been discovered in [|| by the obser- 
vation of the strong pics in the mass spectra of sodium clusters. The enhanced stability 
of some clusters, the so-called magic clusters, was explained by the closure of shells of de- 
localized electrons. A simple physical model describing electronic shell structure of metal 
clusters has been developed within the jellium approximation (see, e.g., ||) by analogy with 
the shell model of atomic nuclei (see, e.g., The jellium model is very successful for the 



simple alkali metals (Na, K), for which one electron per atom is delocalized [Ty, [17|, |18 



The jellium model electronic shell closures for alkali-metal clusters define the magic numbers 
N = 8, 20, 34, 40, 58, 92 that are in a good agreement with experiment. Note that the 
jellium model can be generalized by accounting for the collective ion background vibration 
dynamics fl9| , and be used as a very appropriate framework for the investigating collision 
processes involving atomic clusters [f2l|j . 

Clusters of divalent metals are expected to differ from the jellium model predictions at 
least at small cluster sizes. In this case, bonding between atoms is expected to have some 
features of the van der Waals type of bonding, because the electronic shells in the divalent 
atoms are filled. Thus, clusters of divalent metals are very appropriate for studing non-metal 
to metal transition, testing different theoretical methodologies and conceptual developments 
of atomic cluster physics. However, relatively little work was done so far on the exploration 
of the alkali-earth metal clusters in comparison with that for the alkali-metal clusters; see, 



e.g., [H [22] and references therein. 

Previous theoretical studies of the magnesium cluster properties have been performed 
using pseudopotential approximation for the treatment of inner electrons in a cluster and 
the density-functional theory for the description of outer shell electrons. The electronic 
properties, geometry and stability of small Mg metal clusters with the number of atoms 
N < 7 have been investigated in |23], |24| using the pseudopotential local-spin-density ap- 
proximation. The geometrical structure and bonding nature of Mg^ clusters with N up to 



13 have been studied in |25|] using the density-functional molecular-dynamics method. The 
size evolution of bonding in magnesium clusters Mg^ with iV = 8 — 13, 16, 20 have been 
studied in using the local-density approximation that accounts for gradient corrections. 
Structural and electronic properties of small magnesium clusters (N < 13) were studied in 
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|27j using a first-principles simulation method in conjunction with the density-functional 
theory and the generalized gradient correction approximation for the exchange-correlation 



functional. It was shown [^7j that the metallization in magnesium clusters has a slow and 
nonmonotonic evolution, although, also jellium-type magic clusters were observed [ 55 , |26 



In order to extend such calculations to larger systems, symmetry restricted methods have 
been developed. The spherically-averaged-pseudo-potential scheme with the local and non- 
local pseudopotentials has been used for the investigation of the electronic structure and 
shell closures of spherical Mg N clusters up to N = 46 [p8fl . 

Recently, new experimental data for Mg clusters have been obtained, indicating the most 
intensive pics in the mass spectra at N = 5, 10, 15, 18, 20, 25, 28, 30, 35, 40, 47, 56, 59, 
62, and 69 ||29|| . These numbers deviate from the sequence of magic numbers which were 
obtained for simple alkali metal clusters, and cannot be reproduced within simple jellium 



models. This fact was interpreted in [^, |3(J within the spherical shell model by diving of 
the high angular momentum states down through the states with lower I. 

In the present work we investigate the optimized ionic structure and the electronic prop- 
erties of neutral and singly charged magnesium clusters within the size range N < 21. 
We calculate binding energies per atom, ionization potentials and energy gaps between 
the highest occupied and the lowest unoccupied molecular orbitals. Our calculations are 
based on ab initio theoretical methods invoking the density-functional theory and system- 
atic post-Hartree-Fock many-body theory accounting for all electrons in the system. The 
results obtained are compared with the available experimental data and the results of other 
theoretical works. 

The atomic system of units, |e| = m e — H = 1, has been used throughout the paper, 
unless other units are indicated. 



II. THEORETICAL METHODS 



Our calculations have been performed with the use of the Gaussian 98 software package 



31JI . We have utilized the 6 — 311G(d) basis set of primitive Gaussian functions to expand 



the cluster orbitals [32| . 

The cluster geometries have been determined by finding local minima on the multidi- 
mensional potential energy surface for a cluster. We have taken into account all electrons 
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available in the system, when computing the potential energy surface. With increasing clus- 
ter size, such calculations become computer time demanding. In this work, we limit the 
calculations by the cluster size N = 21. 

The key point of calculations is fixing the starting geometry of the cluster, which could 
converge during the calculation to a local or the global minimum. There is no unique way for 
achieving this goal with Gaussian 98 In our calculations, we have created the starting 
geometries empirically, often assuming certain cluster symmetries. Note, that during the 
optimization process the geometry of the cluster as well as its initial symmetry sometimes 
change dramatically. All the characteristics of clusters, which we have calculated and present 
in next section, are obtained for the clusters with optimized geometry. 

In this work we concentrate on the systematic exploration of the properties of magnesium 
clusters using the density-functional theory based on the hybrid Becke-type three-parameter 
exchange functional |33[] paired with the gradient-corrected Lee, Yang and Parr correlation 
functional (B3LYP) |34], 35fl , as well as the gradient-corrected Perdew-Wang 91 correlation 
functional (B3PW91) |5f| [37]]. The important feature of the density-functional method 
consists in the fact that it takes into account many-electron correlations via the phenomeno- 
logical exchange- correlation potential. However, so far, there has not been found the unique 
potential, universally applicable for different systems and conditions. As a result there 
are many different parameterizations for the exchange-correlation potential valid for special 
cases. 

Alternatively, we use a direct ab initio method for the description of electronic properties 



of metal clusters, which is based on the consistent post-Hartree-Fock many-body theory |38 



In the present work, we apply the M0ller-Plesset perturbation theory of the fourth order 
(MP 4). Based on the fundamental physical principles being free from any phenomenological 
parameters, this model can be refined by extending the quality of the approximations, while 
the physical meaning of the effects included is clearly demonstrated. Thus, often such 
an approach predicts more accurate and reliable characteristics of metal clusters than the 
density-functional theory. 

In the present work we use both different theoretical schemes for calculations taking 
advantage of the clear physical meaning and reliability of the post-Hartree-Fock perturbation 
theory and the numerical efficiency of the density-functional methods. 
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III. NUMERICAL RESULTS AND DISCUSSION 
A. Geometry optimization of Mgjy and Mg^ clusters 




FIG. 1: Optimized geometries of the neutral magnesium clusters Mgi — Mgi\ calculated in the 
B3PW91 approximation. The interatomic distances are given in angstroms. The values in brackets 
correspond to the B3LYP results. The label above each cluster image indicates the point symmetry 
group of the cluster. 

The optimization of the cluster geometries has been performed with the use of the 
B3PW91 and B3LYP methods. For small magnesium clusters with number of atoms 
N < 11, we have also used the ab initio MP A method in addition to density-functional 
calculations. With the growth cluster size the ab initio MP4 calculations become more 
and more computer time demanding, therefore we have not performed such calculations for 
magnesium clusters with the number of atoms iV > 12. The detail comparison of the re- 
sults obtained by the density-functional and ab initio perturbation theory methods as well 




FIG. 2: The same as in Fig.[T] but for singly-charged magnesium clusters Mg 2 — Mg^i- 
as their comparison with the results of other works is given below, see section |lllfj|. This 



comparison allows us to conclude that for magnesium clusters the B3PW91 method is more 
reliable and accurate in comparison with the B3LYP one. 

The results of the cluster geometry optimization for neutral and singly charged magne- 
sium clusters consisting of up to 21 atoms are shown in figures [I] and Q respectively. 

Magnesium clusters possess various isomer forms those number grows dramatically with 
increasing cluster size. In figures p] and 0, we present only the lowest energy configurations 
optimized by the B3PW91 method. The interatomic distances are given in angstroms. The 
values in brackets correspond to the interatomic distances obtained by the B3LYP method. 

Figure [l] shows that the neutral magnesium clusters form the compact structures, max- 
imizing the coordination number. The Mg 2 dimer is weekly bound possessing the binding 
energy per atom 0.039 eV/atom and the bond length 3.609 A, which is in a good agreement 



with the experimental results of Ref. p9j, where the values 0.025 eV/atom for the binding 
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energy and 3.89 A and for the bond length have been reported. The lowest energy state for 
Mg 3 is the equilateral triangle, and for Mg 4 is a regular tetrahedron. As we discuss below, 
the Mg4 cluster is relatively more stable and compact, as compared to the neighbouring 
clusters. The Mg§ cluster has a structure of slightly elongated triangular bipyramid, while 
Mg§ consists of three pyramids connected by their faces, Mg-j is a pentagonal bipyramid, 
and Mg s is a capped pentagonal bipyramid. These geometrical structures are in a good 
agreement with the results of Ref. [^5| . 



It is worth to note that the optimized geometry structures for small neutral magnesium 



clusters differ significantly from those obtained for sodium clusters (see, e.g., p2|, [40], [41| 
and references therein). Thus, the optimized sodium clusters with N < 6 have the plane 
structure. For Na^, both plane and spatial isomers with very close total energies exist. The 
planar behavior of small sodium clusters has been explained as a result of the successive 



filling of the la and In symmetry orbitals by delocalized valence electrons |4(J, which is 
fully consistent with the deformed jellium model calculations [18|]. Contrary to the small 
sodium clusters, the magnesium clusters are tri-dimensional already at N = 4, forming the 
structures nearly the same as the van der Waals bonded clusters. 

Starting from Mgg a new element appears in the magnesium cluster structures. This is 
the six atom trigonal prism core, which is marked out in figure [|. The formation of the 
trigonal prism plays the important role in the magnesium cluster growth process. Adding 
an atom to one of the triangular faces of the trigonal prism of the Mgg cluster results in the 
Mgio structure, while adding an atom to the remaining triangular face of the prism within 
the MgiQ cluster leads to the structure of Mgu, as shown in figure |l|. 

Further growth of the magnesium clusters for 12 < iV < 14 leads to the formation of the 
low symmetry ground state cluster. In spite of their low symmetry, all these clusters have 
the triagonal prism core. The structural rearrangement occurs for the Mg 15 cluster, which 
results in the high symmetry structure of the two connected Mgg clusters. 

Starting from Mg\§ another motif based on the hexagonal ring structure which is marked 
out in figure |] dominates the cluster growth. Such a ring is the basic element of the hexagonal 
closest-packing (hep) lattice, as one can see in figure |3], in which the primitive cell for the 
magnesium hep lattice is presented. Thus, N = 15 is the turning-point in the formation of 
the hep lattice for magnesium. 

Vectors a, b and c in figure ^| show the primitive cell axes of the hep lattice. For bulk 
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FIG. 3: Primitive cell for magnesium hep lattice. For bulk magnesium a = b = 3.21 A and c = 5.21 

A@. 

magnesium a = b = 3.21 A and c = 5.21 A El. The fundamental characteristic for the 
hexagonal closest-packing of spheres is the value of ratio c/ a, which is equal to a/8/3 ~ 1.633 
for ideal hep lattice. The bulk magnesium with c/a = 1.62 is very close to ideal hep structure 

The distinct three- layered structure of Mg-clusters with N > 18 based on the hexagonal 
rings allows one to determine the averaged values of the primitive axes (c) and (a). Table | 
demonstrates that the calculated values (c) and (a) and their ratio for magnesium clusters 
with N > 18 are very close to the corresponding values for bulk magnesium. 

Figure |2] shows the optimized geometries of singly- charged cationic magnesium clusters. 
The ground state geometries of the cationic magnesium clusters are not very different from 
those obtained for the neutral parent clusters with the exception of Mg^ and Mg£, the 
equilibrium geometries of which are linear chains. Below, we discuss the stability of the 
linear chain isomers for the magnesium clusters (neutral and singly-charged) within the size 
range considered. 

In figure [|, we present the average bonding distance, (d), calculated within the B3PW91 
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TABLE I: The averaged values of the primitive axes and its ratio for the hep lattice element 

for magnesium clusters with N > 18 calculated within the B3PW91 approximation. Values in 
brackets correspond to singly-charged magnesium clusters. 

Mg 18 Mg 19 Mg 20 Mg 21 Mg bulk, @|] 

(c), A 5.08 (5.42) 5.47 (5.37) 5.48 (5.23) 5.56 (5.23) 5.21 

(a), A 3.14 (3.19) 3.05 (3.22) 3.20 (3.21) 3.20 (3.21) 3.21 

(c)/(a) 1.62 (1.70) 1.79 (1.67) 1.71 (1.63) 1.74 (1.63) 1.62 



approximation for neutral and singly-charged magnesium clusters. When calculating the 
average bonding distance in a cluster, interatomic distances smaller than 4.1 A have only 
been taken into account. The bulk limit for the magnesium hep lattice indicated in 
figure by horizontal dashed line. 




FIG. 4: The average bonding distance as a function of cluster size for neutral and singly-charged 
magnesium clusters. Stars present the results of the work by Akola et al [^?]]. The horizontal 
dashed line indicates the bulk limit for the hep lattice [^] . 
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Figure ^ shows how the average bonding distance evolve with increasing cluster size. It 
is clearly seen that the dependence of the average bonding distance on cluster size has es- 
sentially non- monotonous oscillatory behavior. For Mg 2 , the bonding distance calculated 
within the B3PW91 method is equal to 3.609 A, which is in a good agreement with the 
experimental result 3.891 A of Ref. P5fl . The appearance of the minima in the size de- 
pendence of the average bonding distance shows that Mg 4 , Mg w , and Mg 2 o clusters (8, 20 
and 40 valence electrons, respectively) are more tightly packed than their neighbours. This 
behavior can be interpreted by the influence of electronic shell effects on the geometrical 



structure of magnesium clusters. It supports the conclusion of Ref. ||4j that electronic shell 
effects can enhance the stability of geometric structures resulting from dense ionic packing. 

Additional minimum in the dependence of the average bonding distance on N arrises at 
N = 15. At this iV a considerable rearrangement of the cluster geometry take place as it is 
seen from figure |]. Indeed, starting from the M g 15 cluster the three-layered structure based 
on the hexagonal ring is formed. It is clearly seen in figure [| that for iV > 15 the average 
bonding distance for magnesium clusters approaches the bulk limit. 

The evolution of the average bonding distance with cluster size differs for magnesium 
clusters from that for sodium. For neutral sodium clusters, one can see odd-even oscillations 
of (d) atop its systematic growth and approaching the bulk limit ||22|| . These features have 
the quantum origin and arise due to the spin coupling of the delocalized valence electrons. 
For magnesium clusters, the average bonding distance depends on size non-monotonically, 
with minima for the Mg^, Mgio, Mgi%, and Mg 2 o clusters. Such an irregular behavior is 
induced by both the closure of electronic shells of the delocalized electrons and structural 
rearrangement s . 

Manifestation of the magic numbers in the dependence of the average bonding distance 
on cluster size coinciding with the spherical jellium model magic numbers does not imply, 
however, the rapid metallization of magnesium clusters. To investigate the transition of van 
der Waals to metal bonding in magnesium clusters it is necessary to explore in detail the 
evolution of their electronic properties. Below we perform such analysis in detail. 

Dashed line in figure |] shows the average bonding distance as a function of cluster size 
calculated for singly-charged magnesium clusters. Figure f| demonstrates the essential dif- 
ference in the cluster size dependence of (d) for the cationic and neutral magnesium clusters 
with N < 6. The small cationic magnesium clusters are more compact in comparison with 
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the corresponding neutral clusters. For example, for Mg£ the bonding distance is equal to 
3.044 A, which is much less than in the case of Mg 2 . This phenomenon has a simple physi- 
cal explanation: the removed electron is taken from the antibonding orbital. The fact that 
cationic magnesium clusters are more stable than the parent neutral and anionic clusters 



has been already noted in ]2j] 



Within the size range N >7, the average bonding distances for single- charged and neutral 
magnesium clusters behave similarly. The absolute value of (d) for single-charged clusters 
is slightly larger in this region of N. 

Figure f| demonstrates the good agreement of our results with the dependence of (d) on 
N calculated in |27|] for neutral Mg-clusters within the size range iV < 13. 

B. Binding energy per atom for Mgjy and Mg^ clusters. 

The binding energy per atom for small neutral and singly-charged magnesium clusters is 
defined as follows: 

E b /N = E x - E N /N (1) 
E+/N = ((N - 1)^ + E+ - E+) /N, (2) 

where En and E N are the energies of a neutral and singly-charged N-particle atomic cluster, 
respectively. E± and Ef are the energies of a single magnesium atom and an ion. 

Figures |^ and |6] show the dependence of the binding energy per atom for neutral and 
singly-charged clusters as a function of cluster size. The energies of clusters have been 
obtained using the B3LYP, B3PW91 and MP4 methods. Calculations of the binding en- 
ergies have been performed by different theoretical methods and with the use of different 
exchange-correlation functionals for the sake of comparison of their accuracy and computa- 
tion efficiency. In figure |5] filled rhombus, crossed rhombus and opened pentagons show the 
result of calculations by Kumar et al [^J, Reuse et al [^4]] and Delaly et al |26j respectively. 



These calculations have been performed within the Hohenberg-Kohn-Sham local-density ap- 
proximation using the Perdew and Zunger [45[] parameterization of the Ceperley and Alder 



46| data for the exchange correlations. Crossed circles and stars present the results of Delaly 



et al [^{J and Akola et al [^7]] derived with the use of the gradient-corrected approximation 
47, 4g| and the PBE parameterization of the gradient-corrected exchange-correlation energy 
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functional E9| respectively. 




FIG. 5: Binding energy per atom for neutral magnesium clusters as a function of cluster size. 
Squares, circles and lower triangles represent the binding energies per atom calculated by the 
B3LYP, B3PW91 and MPA methods respectively. Stars, rilled rhombus and crossed rhombus 



show the results of the works by Akola et al [27], Kumar et al [25], and Reuse et al [24] respectively. 



Opened pentagons and crossed circles show the result of Delaly et al [26] obtained with the use 
of the LDA and gradient-corrected (GC) LDA methods respectively. Labels indicate the point 



symmetry group of the isomers represented. Their geometries one can find in section III A 



Figure |5] shows that, although, the qualitative behavior of the binding energy per atom 
calculated within different approaches is similar, the quantitative discrepancy between the 
curves is rather considerable. This is a result of different accounting for the gradient cor- 
rections to the local-density exchange correlation interaction within different methods. The 
gradient corrections have been shown to provide a systematic improvement in the computed 



properties of magnesium clusters [p6| . The difference in the binding energy per atom for 
neutral magnesium clusters with AT < 21 calculated with the use of the gradient corrected 
B3LYP and B3PW91 methods reaches 0.35 eV. The reason for this difference is in the 
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FIG. 6: The same as in Figj5] but for singly-charged magnesium clusters. 

different way of accounting for many-electron correlations within the B3LYP and B3PW91 
methods. To explore what type of parameterization of the exchange-correlation energy is 
more reliable for magnesium clusters we have used the post-Hartree-Fock M0ller-Plesset 
perturbation theory. This method is free of phenomenological parameters and can be used 
as a criterion for checking the reliability of various density-functional theory schemes. The 
disadvantage of the perturbation theory approach consists in the fact that it leads to the 
dramatic growth of the computational costs with increasing the number of electrons in the 
system in comparison with that for the density-functional theory calculations. Therefore, 
we have used the MP A method only for clusters with the number of atoms N < 11. 

Figure |5| shows that the results of the MP4 theory are in a reasonable agreement with 
those derived by the B3PW91 method. This comparison demonstrates that for magnesium 
clusters simulations the B3PW91 method is more reliable than the B3LYP one. Our results 
derived within the B3PW91 and MP4 approximations are in a good agreement with those 
from Ref. |[Z§ £7]]. 



We now discuss the behavior of the binding energy as a function of cluster size for both 



14 




FIG. 7: Second differences of total energy for neutral, A 2 E N = E N+1 — 2E N + E N _ 1 , and singly- 



charged, A 2 E+ = E+ - 2E+ + E+_ v magnesium clusters. Stars show the result of the work by 



Akola et al [g7|, filled rhombus by Kumar et al |2{|], and crossed rhombus by Reuse et al [24 |. 



neutral and singly-charged magnesium clusters. For neutral magnesium clusters, the binding 
energy per atom increases steadily with the growth cluster size. The local maxima of E^/N 
at iV = 4, 10 and 20 correspond to the most stable configurations of the magic magnesium 
clusters possessing N e i = 8, 20 and 40 valence electrons respectively. This behavior is in 
agreement with the simple spherical jellium model. The analysis of the second differences 
of the binding energy (see Fig.[7|) confirms this conclusion and makes a hint about relative 
stability of the Mg?, Mgi 3 , Mgi 5 , and Mgu clusters, in addition to the magic clusters Mg 4 , 
Mg w and Mg 2 o . The additional magic numbers can be explained within the deformed 



jellium model accounting for spheroidal deformations of the cluster core (see, e.g., |16|, [17], [18 
and references therein). For a spheroidal jellium cluster, the orbital angular momentum does 
not remain a good quantum number characterizing the valence electrons energy levels. In 
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this case, the energy levels are characterized by the projection of the angular momentum A 
on the principal axis and by the parity of the wave function. Thus, the energy levels with 
A = are twofold degenerated on the projection of the electron's spin, while those with 
A 7^ are fourfold degenerated both on the projection of the electron spin and on the sign of 
the projection A on the principal cluster axis. The deformed jellium clusters having closed 
electronic subshells possess the enhanced stability. Therefore, in addition to the spherical 
magic clusters with 8, 20, 40 etc valence electrons, the deformed jellium clusters with 6, 10, 
14, 18, 22, 26, 30, 34 etc. valence electrons turn out to be relatively stable. This fact leads 
to the following additional magic numbers 3, 5, 7, 9, 11, 13, 15, 17 for the jellium magnesium 
clusters. Some of these numbers, such as 3, 5, 9, 11 precede or follow the spherical magic 
numbers 4, 10, 20 and as a result of that become masked and are not that pronounced in 
the second differences analysis. 

For singly- charged magnesium clusters, the binding energy per atom as a function of 
cluster size is essentially non-monotonous. The local maxima of the binding energy for 
the Mg%, Mg£, Mg± , Mgf 2 , Mgf 5 and MgJ, clusters indicate their enhanced stability. 
Figure [1] shows second differences of the total energy for singly-charged magnesium clusters. 
This figure demonstrates the enhanced stability of the mentioned cluster ions and the Mgf 7 
cluster. 

The sequence of magic numbers for singly charged magnesium clusters differs from that 
for neutral clusters. This happens because singly charged magnesium clusters always possess 
odd number of valence electrons and, thus, always contain open electronic shells. For neutral 
magnesium clusters, situations of both close and open electronic shells are possible. The 
enhanced stability of a Mg-cluster ion arises, when its electronic configuration has one 
hole in or an extra electron above the filled shells. Thus, the cluster ions Mg£, Mgf ± 
and Mg%i contain one extra electron over the complited spherical electronic shells, while 
the clusters Mg\ , Mgf and Mg^ have a hole in the spherical outer electronic shell. Our 
results presented in figures |6] and [7| demonstrate that the cluster ions Mg£, Mgf and MgJ, 
turn out to be more stable than their neighbors. We note that the alteration of the magic 
number from N = 4 for neutral Mg-clusters to N = 5 for Mg-cluster ions happens because 
the electronic configuration containing an extra electron becomes more favorable for Mg£. 
This is not the case for the Mgf and Mg 20 clusters, those outer electronic configurations 
contain a hole. 
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The Mg-cluster mass spectra have been recorded in [29[ indicating the enhanced stability 
of the clusters with N = 5, 10, 15, 18 and 20. In that work the role of the cluster ionization 



was not reliably clarified |29| and thus the charge state of the clusters was not reliably 
determined. As a result, the observed magic numbers sequence should be a combination of 
the magic numbers sequences for neutral and singly-charged cluster ions. Thus, N = 5 is 
the ionic magic number, N = 10, 15 and 20 are the magic numbers manifesting themselves 
clearly for both neutral Mg-clusters and Mg-cluster ions. The second differences are positive 
and relatively large for N = 13 (neutral clusters) and N = 12 (singly-charged cluster ions). 
Possibly, the interplay between neutral clusters and ions make these numbers masked in 
experiment. The second differences are also positive for N = 7 for neutral Mg-clusters and 
for iV = 3 for Mg-cluster ions, although the enhancement for those numbers have not been 
experimentally observed. We explain this fact by possible suppression of the experimental 
signal in the region of small N and relatively small values of the second differences in the 
mentioned cases. 

The potential energy surface for a cluster becomes more and more complicated with 
increasing cluster size. The magnesium clusters are not an exception. Figure § demonstrates 
this fact where we present the binding energies per atom calculated for a variety of isomers 
of neutral magnesium clusters. The corresponding point symmetry groups and the accurate 
values of the total energies calculated within the B3LYP and B3PW91 approximations are 
presented in Appendix in tables [TT] and f[TT] respectively. Most of the isomer configurations 
have been obtained using the B3LYP method, while the B3PW91 method has been used 
for the exploration of the ground state energy isomers, as well as for the linear and ring-like 
isomer structures. 

Squares in figure |8] correspond to the most stable clusters possessing the minimal total 
energy. Among the variety of isomers, presented in figure ||, we mark certain groups of 
isomers with the fixed symmetry. So, circles present the linear chains point symmetry 

group) and the upper triangles correspond to the rings of N atoms (D^h point symmetry 
group). It is an interesting fact that among the multitude of the isomers of neutral magne- 
sium clusters the linear chains and rings are always stable. We pay a particular attention 
to these structures because of their possible applications in nano-technology. Extracting 
these isomers and putting them on a substrate one can produce one-atom wide quantum 
wires. The linear chains and rings of atoms are also very interesting from the theoreti- 
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FIG. 8: Binding energy per atom for a variety of isomers of neutral magnesium clusters as a 
function of cluster size. The corresponding point symmetry groups and the accurate values of the 
total energies are presented in Appendix in tables [ll] and III. Numbers near some lines show the 
number of found isomers with the corresponding close energies. 
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cal point of view, because with their help one can investigate the transition from one- to 
two-dimensional structures. For linear chains, the binding energy per atom increases slowly 
with the growth the number of atoms, while in the case of rings the value of E^/N has the 
prominent odd-even oscillatory behavior. 

This behavior arises as a result of successive filling the a- and 7r-symmetry orbitals by 
valence electrons in magnesium linear chains and rings. Indeed, because of its symmetry 
the one-dimensional linear chain of N magnesium atoms has the following configuration of 
valence electrons: la 2 , 2a 2 , 3a 2 , Na 2 . Therefore for any N it has the closed electronic 
shell structure. This fact explains the monotonous growth with N of the linear chain binging 
energy, and its relative saturation in the region iV > 10. 

The molecular orbitals for the structure of the ring-type have to be aligned with the plane 
of the ring. Such orbitals are fourfold degenerated due to symmetry reasons. The Mg 2 dimer 
has four valence electrons that occupy spherically-shaped la 2 and prolate-like 2a 2 orbitals. 
The Mg% trimer has six valence electrons, two of them occupy la 2 state, while the remaining 
four electrons fill the fourfold degenerated l7r 4 orbital aligned with the plane of the trimer. 
With increasing the number of magnesium atoms in the ring, the valence electrons continue 
to occupy fourfold degenerated orbitals aligned with the plane of the ring. Therefore, in the 
magnesium ring-like isomers with the odd number of atoms all electronic shells are closed, 
while the isomers with the even number of atoms have the open electronic shell. This fact 
results in the enhanced stability of the magnesium rings with an odd number of atoms 
N = 3, 5, 7, 9, ... etc. and explains the odd-even oscillatory behavior of the binding energy 
for the magnesium rings. 



C. Ionization potentials and HOMO-LUMO gaps 

Let us now consider how the ionization potential of magnesium clusters evolves with 
increasing cluster size. The ionization potential of a cluster V{ is equal to the difference 
between the energies of the corresponding cluster ion and the neutral cluster, Vi = E^ — E N . 
Figure ||| shows the dependence of the adiabatic V°~ diab (i.e. the geometry relaxation of the 
ionized cluster is taken into account) and vertical V^ ert (i.e. the cluster geometry is frozen 
during the ionization process) ionization potential on N. We compare our results derived 



by the B3PW91 method with theoretical data from Ref. [E7fl and [24| and with the bulk 
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limit, V„ bulk = 3.64 eV, taken from M 



Both the vertical and adiabatic ionization potentials evolve non-monotonously with in- 
creasing cluster size. Figure ^| shows that ionization potential of magnesium clusters steadily 
but rather slow decreases towards the bulk limit. This evolution is neither rapid nor 
monotonous process. In order to exclude the influence of the cluster geometry rearrange- 
ment, we first consider the vertical ionization potential. The size dependence of the vertical 
ionization potential has a prominent maximum at N = 4 followed by a sharp decrease. Such 
a behavior of the ionization potential is typical for the jellium model, predicting maxima 
in the size dependence of the ionization potential at the magic numbers corresponding to 
the clusters with closed electronic shells. Our data are in a good agreement with the results 



of Ref. 1241 , but contradict to those reported in Ref. [^7| for the Mg% and Mg± clusters. 



In [27] the appearance of the deep minimum in vertical ionization potential at N = 4 was 
explained as a result of a stronger charge derealization in the Mg 4 cluster in comparison 
with its neibours. 

We note that the peculiarities in the ionization potential dependence on N correlate 
with the magic numbers that appear for the singly-charged magnesium clusters. Indeed, 
the minima in V^ ert correspond to the maxima in E^/N for Mg-cluster ions (see Fig. ^j. 
This fact has a simple explanation. The ionization potential of a cluster is equal to the 
difference between the energies of the corresponding cluster ion and the neutral cluster. 
For neutral Mg-clusters, the binding energy as a function of N increases steadily with the 
growth cluster size, while for Mg-cluster ions - irregularly. Thus, their difference mimics 
all the irregularities that appear in the binding energy dependence on N for singly-charged 
magnesium clusters. 

For N > 6, the vertical ionization potential changes slowly with increasing cluster size. 
This process is characterized by the irregularities that originate due to the influence of the 
cluster geometry on the jellium-type electronic structure of Mg-clusters. 

Indeed, the shape of a jellium cluster is defined by its electronic structure. Thus, the 
closed shell jellium clusters are spherical, while clusters with opened electronic shells are de- 
formed due to the Jahn- Teller distortions. The jellium picture works fairly well for sodium 
clusters. The ionization potential of sodium clusters drops rapidly and systematically at the 
electronic shell closures. The iV-dependence of the ionization potential has prominent, reg- 
ular odd-even oscillations (see, e.g., |TR and references therein). Magnesium clusters are 
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FIG. 9: Adiabatic, V i adtab , and vertical, V^ ert , ionization potential for small magnezium clusters. 
Stars and crossed rhombus show the result of the work by Akola et al [27] and by Reuse et al [24] 
respectively. 
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different. As discussed in section [III A| , the evolution of the Mg-cluster geometry is closely 
connected with the formation of elements of the hep lattice cell. Although, the electronic 
shell effects clearly manifest themselves in the formation of the Mg-cluster geometry, they 
do not determine it completely. Another words, there is an interplay of the jellium and the 
hep lattice factors in the formation of the Mg-cluster geometry and the electronic properties 
such as the binding energy and the ionization potential. 

The adiabatic ionization potential dependence that is shown in figure |9] exhibits quali- 
tatively the same behavior as the vertical one, however, has more pronounced irregularities 
due to geometry rearrangements of the ionized clusters. 



Figure [10] shows the gap E g between the highest occupied and the lowest unoccupied 
molecular orbitals (HOMO-LUMO gap) for the Mg-clusters as a function of cluster size. 
For the sake of comparison, we have also calculated the HOMO-LUMO gap for the sodium 
clusters and present it in figure |H]. Calculations have been performed using the B3PW91 
and B3PLYP methods. The geometries of neutral sodium clusters have been taken from 
|22fl . For the small magnesium clusters with iV < 13, we compare our results with those 



presented in Ref. J27] . 

The size dependence of E g for neutral sodium clusters has an oscillatory behavior with 
local maxima at iV = 6, 8, 10, 14 and 20. These maxima correspond to the electronic shell 
closures in a full accordance with the deformed jellium model. The local maximum in the 
size dependence of E g at N = 12 and the shift of the local maximum from N = 18 to N = 17 
are the consequences of triaxial deformations ||22|| . Thus, the triaxial deformation leads to 
the splitting of the fourfold degenerated highest occupied orbital on two twofold degenerated 
orbitals. As a result of that the additional shell closure at the N e i = 12 appears. 

For Mg-clusters, the evolution of the HOMO-LUMO gap with the growth cluster size 
differs from that for iVa-clusters. The gap E g calculated for magnesium clusters shows the 
oscillatory behavior accompanied by the gradual decrease in the absolute value. Maxima 
in this dependence at N = 4, 10 and 20 correspond to the magic numbers of the spherical 
jellium model (N e i = 8, 20 and 40 respectively). The similar feature also does exist for 
iVa-clusters at iV = 8 and 20. Additional variation of E g appears both due the subshell 
closures and the cluster structural rearrangements. 

We note that the HOMO-LUMO gap remains rather large even for the clusters with 
N > 15 possessing elements of the hep lattice of the bulk Mg. This fact confirms the 
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FIG. 10: Gap between the highest occupied and the lowest unoccupied eigenstates for the Mg and 
Na clusters as a function of cluster size. Circles and squares represent the HOMO-LUMO gap 
calculated by the B3PW91 and the B3LYP methods respectively. Stars show the result of the 
work by Akola et al . 
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conclusion on the slow and non-monotonous evolution of metallic properties in Mg clusters. 

IV. CONCLUSION 

The optimized geometries and electronic properties of neutral and singly- charged mag- 
nesium clusters consisting of up to 21 atoms have been investigated using the B3PW91, 
B3LYP and MP4 methods accounting for all electron in the system. The detailed com- 
parison of the results of the phenomenological B3PW91 and B3LYP density-functional 
methods with the results of the systematic ab initio post-Hartree-Fock many-body theory 
leads us to the conclusion that the B3PW91 method is more reliable for Mg-cluster simu- 
lations than the B3LYP one. 

We have investigated the size evolution of the Mg-clusters geometry. It has been shown 
that starting from Mgi§ the hexagonal ring structure determines the cluster growth, which 
is the basic element of the hep lattice for the bulk magnesium. 

We have investigated the electronic properties of magnesium clusters. It has been shown 
that the electronic shell effects and jellium-like behavior clearly manifest themselves in the 
formation of geometrical properties, however, the shell effects do not determine the geometry 
of Mg clusters completely. We have demonstrated that due to the interplay of the jellium and 
the hep lattice factors the electronic properties of magnesium clusters possess irregularities 
which can not be explained within the simple jellium model. It has been shown that the 
metallic evolution of magnesium clusters is slow and non-monotonous process. 

The results of this work can be extended in various directions. One can use the similar 
methods to study structure and properties of various types of clusters. It is interesting to 
extend calculations towards larger cluster sizes and to perform more advanced comparison 
of model and ab initio approaches. A lot of novel problems arise, when considering collisions 
and electron excitations in the clusters with the optimized geometries |2"I |. These and many 



more other problems on atomic cluster physics can be tackled with the use of methods 
considered in our work. 
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TABLE II: Total energies and the point symmetry groups for a variety of isomers of neutral 



magnesium clusters. Calculations have been done by the B3LYP method. 



N 


Point group 


Energy (a.u.) 


N 


Point group 


Energy (a.u.) 


N 


Point group 


Energy (a.u.) 


1 




-200.0931 


9 


Dooh 


-1800.8472 


16 


d.Cs 


-3201.6766 


2 


D^h 


-400.1868 


10 


c 3v 


-2001.0336 




C.Cs 


-3201.6713 


3 


D 3h 


-600.2840 




Civ 


-2001.0239 




b.Ci 


-3201.6673 




Doch 


-600.2807 




T d 


-2001.0225 




a.Ci 


-3201.6646 


4 


T d 


-800.3938 




Ci 


-2001.0120 




T d 


-3201.6594 




D 2h 


-800.3792 




D 3h 


-2001.0114 




a.Cs 


-3201.6579 




D^n 


-800.3750 




C 2v 


-2001.0098 




b.Cs 


-3201.6572 




D 4h 


-800.3748 




D 4d 


-2001.0078 




Dooh 


-3201.5082 


5 


D 3h 


-1000.4907 




Dooh 


-2000.9417 


17 


D id 


-3401.7953 




C 1 4 v 


-1000.4787 




DlQh 


-2000.9322 




Cs 


-3401.7861 




Td 


-1000.4768 


11 


D 3h 


-2201.1348 




D 3h 


-3401.7052 




C'2v 


-1000.4766 




D llh 


-2201.0430 




Dooh 


-3401.6026 




D 5h 


-1000.4751 




Dooh 


-2201.0362 


18 


b.C 2v 


-3601.9095 




D 2D 


-1000.4743 


12 


a.Cs 


-2401.2366 




C.Cs 


-3601.9026 




D 2h 


-1000.4738 




b.Cs 


-2401.2304 




c 2 


-3601.8764 






-1000.4694 




c 2v 


-2401.2178 




b.Cs 


-3601.8752 


6 


C 2v 


-1200.5898 




D 6h 


-2401.1313 




a.Cs 


-3601.8716 




D 2h 


-1200.5897 




Dooh 


-2401.1308 




a.C 2v 


-3601.8711 






-1200.5851 




D 12h 


-2401.1278 




D 5h 


-3601.8608 






-1200.5815 


13 


b.Ci 


-2601.3438 




D 6h 


-3601.7022 




D 3h 


-1200.5765 




Cs 


-2601.3308 




Dooh 


-3601.6970 




o h 


-1200.5720 




a.Ci 


-2601.3302 


19 


c 2v 


-3802.0292 






-1200.5639 




c 3v 


-2601.3221 




c 3v 


-3801.9867 






-1200.5638 




Ih 


-2601.2904 




D 5h 


-3801.9706 


7 


D 5h 


-1400.6933 




D 6h 


-2601.2401 




D 6h 


-3801.8449 




c 3 


-1400.6908 




Dooh 


-2601.2253 




Dooh 


-3801.7925 




C 3v 


-1400.6854 


14 


Ci 


-2801.4485 


20 


Cs 


-4002.1503 




D 7h 


-1400.6645 




b.C 3v 


-2801.4397 




Ci 


-4002.1444 




Dooh 


-1400.6583 




C s 


-2801.4407 




d.C 2v 


-4002.1392 


8 


a.C 2v 


-1600.7999 




c 2v 


-2801.4328 




c.C 2v 


-4002.1263 




Cs 


-1600.7976 




o h 


-2801.4115 




b.C 2v 


-4002.1212 




b.C 2v 


-1600.7948 




a.C 3v 


-2801.4134 




T d 


-4002.1026 




T d 


-1600.7854 




Dooh 


-2801.3194 




a.C 2v 


-4002.0990 






-1600.7527 


15 


D 3h 


-3001.5692 




D 6h 


-4001.9764 


9 


D 3h 


-1800.9162 




Cs 


-3001.5627 




Dooh 


-4001.8870 




c 3v 


-1800.9064 




D 6h 


-3001.4594 


21 


Ci 


-4202.2460 






-1800.8540 




Dooh 


-3001.4138 




c 2v 


-4202.2255 
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TABLE III: The same as 


table P but for the P3PW91 method. 


TVT 
IN 


Point group 


Energy (a.u.) 


N 


Point group 


Energy (a.u.) 


1 




-200.0379 


14 


Ci 


-2800.8170 


2 




-400.0788 




flood 


-2800.6007 


3 


D 3h 


-600.1281 




Dl4h 


-2800.5997 






-600.1207 


15 


D 3h 


-3000.9045 


■1 


Td 


-800.1962 




C s 


-3000.8914 






-800.1638 




D 15h 


-3000.6786 






-800.1633 




Dooh 


-3000.6453 


5 


D 3h 


-1000.2450 


16 


d.Cs 


-3200.9666 






-1000.2193 




C.Cs 


-3200.9560 




Dash 


-1000.2063 






-3200.6902 


6 


C 2v 


-1200.2970 




Dl6h 


-3200.6880 




D 2 h 


-1200.2966 


17 


D 4d 


-3401.0612 




C& v 


-1200.2865 




C 3 


-3401.0428 




D 6h 


-1200.2495 




D 17h 


-3400.7719 






-1200.2495 






-3400.7353 


7 


Dsh 


-1400.3604 


18 


b.C 2v 


-3601.1279 




D 7h 


-1400.3091 




C.Cs 


-3601.1221 






-1400.2928 




b.C s 


-3601.0756 


8 


C s 


-1600.4205 




Dcah 


-3600.7806 




a.C2v 


-1600.4193 




D ish 


-3600.7766 




T d 


-1600.4005 


19 


C2v 


-3801.2093 






-1600.3364 




Dim 


-3800.8653 




D 8h 


-1600.3335 






-3800.8260 


9 


D 3h 


-1800.5018 


20 


C s 


-4001.2920 




Dg h 


-1800.4000 




Ci 


-4001.2801 






-1800.3800 




b.C 2v 


-4001.2539 


10 


C 3v 


-2000.5786 




D 2 oh 


-4000.9111 




Td 


-2000.5590 






-4000.8716 






-2000.4367 


21 


Ci 


-4201.3398 






-2000.4238 




C 2 v 


-4201.3201 


11 


D 3h 
D llh 

Doch 


-2200.6322 
-2200.4923 
-2200.4678 








12 


a.Cs 
D 12h 


-2400.6891 
-2400.5321 
-2400.5119 








13 


b.d 
D 13h 


-2600.7561 
-2600.5853 
-2600.5562 
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